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ABSTRACT 

Test runs of the transonic axial compressor test rig at the Naval Postgraduate School, 
Turbopropulsion Laboratory, were conducted in preparation for the installation of a new 
Stage design. Modifications in the cooling air supply to the high speed bearings, and to the 
design of the torque measuring system were completed during subsequent overhaul. 

A case study of the design of the new transonic stage was initiated. This consisted of 
a review of the procedure used in the design, as well as a design comparison. The 
comparison examined the differences between the blades designed for the new stage, which 
was primarily accomplished using a full, three-dimensional, Computational Fluid Dynamics 
code, and blades designed using two-dimensional streamline curvature methods. The axi- 
Symmetric through-flow code, used in the design case study was modified to run on 
workstations at the Naval Postgraduate School, Department of Aeronautics and Astronautics, 


providing students and faculty with a design tool for single or multiple stage axial flow 


compressors. 
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LIST OF SYMBOLS 


Specific Head Rise Parameter 
Specific Heat at Constant Pressure 
Inlet Total Temperature 

Pressure Ratio, Total-to-Total 


Ratio of Specific Heats 


Rotor Wheel Speed at Tip 


I. INTRODUCTION 


The Naval Postgraduate School (NPS) transonic axial compressor design was 
completed by Dr. M.H. Vavra in July of 1968. All of the calculations and drawings 
required for the complete aerodynamic and mechanical design were done by Dr. Vavra 
alone, by hand. From 1968 until the mid 1980’s the compressor served well as a test 
vehicle for the development of innovative flow measurement instrumentation as well as a 
means for graduate students to gain experience in the operation, test, and analysis of high 
speed compressors [Ref. 1-5]. The initial stage design was unique but no real technology 
advance was attained. This was due to two design inaccuracies, discovered after some 
period of testing. In addition to an incorrect assumption of constant through-flow velocity, 
an error was made in the calculations associated with the rotor blade setting angle, and 
this was subsequently built into the blading [Ref. 6]. Consideration was given to twisting 
the fabricated blades to correspond with the shape intended in the design. This was not 
attempted since the design itself had resulted from an inaccurate through-flow prediction. 
Following the completion of a study by Neuhoff [Ref. 7], to measure rotor losses, and to 
separately identify shock and viscous components, the compressor was not operated again 
until the present study was initiated. 

A transonic axial compressor stage designed specifically for the existing NPS 
transonic test rig, shown in Figure 1, was completed at NASA Lewis in 1994, by Nelson 
L. Sanger [Ref. 8]. This new design was sought to perform the functions originally 
intended for the Vavra design, namely, to serve as a test vehicle for instructional purposes, 
and to provide a tool for meaningful research. However, the timing of the project, after a 
period of years in which advanced Computational Fluid Dynamics (CFD) methods for 
turbomachinery had evolved, provided a unique opportunity. The design would provide a 
test of a wholly CFD design approach. The experimental evaluation would provide both a 


validation of the design approach, and an experimental test case for CFD analysis codes. 


The manufactured blading is shown in Figure 2. 
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Figure 1. Transonic Compressor Test Rig. 
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Figure 2. Manufactured Blading; Rotor (Top), Stator (Bottom). 


This report is intended to document the effort to bring the NPS transonic 
compressor test rig back on line, in order to evaluate the newly designed axial stage, and 
to review, for the record, the procedure used in the stage design. The report has two main 
sections. The first section contains a review of the process used in the design of the new 
stage. This section also contains a comparison of the final blade design, resulting from 
CFD analysis, with a blade design using conventional streamline curvature and 
axisymmetric through-flow calculations. The second section is an account of the process 
that was conducted to return the test rig to working order. Appendix A contains a 
description of, AXIDES, a 2-D blade design code, written by Jim Crouse [Ref. 9], which 
was used in the preliminary design and in the production of final fabrication coordinates of 
the new stage. This code has been modified for use at the NPS, Department of 
“Aeronautics and Astronautics, as a tool for advanced axial compressor design. Appendix 


C contains information regarding the test rig disassembly. and reassembly, and the 


operating procedure. 


II. NEW STAGE DESIGN 


A. PROCEDURE 


1. Design Intent 


The object of the design was to produce a state-of-the-art transonic compressor 
Stage, representative of a typical inlet stage, or fan stage. At the time the design was being 
contemplated, the use of CFD codes as analysis tools was increasing. However there was 
little work done using these codes for design. Given that the stage was to be used for 
research and was not intended as an actual aircraft component, the decision was made to 
use CFD as the principle tool in the aerodynamic design. Additional computer codes were 
used in the design including the AXIDES code, however these codes were used mainly in 
supporting roles. The complete package, including the design, testing and measurement of 
the stage constitutes, in the words of the designer, “the ultimate CFD validation 


experiment”. 


2. Preliminary Calculations 


Although the design was intended to be derived from State-of-the-art technology, 
the decision to utilize the existing NPS Turbopropulsion Laboratory (TPL) transonic 
compressor test rig imposed several constraints that were unavoidable. The decision to use 
the existing test rig mechanical design was economic in nature. Retooling the rig to 
produce a desired flow path was determined to be unjustified in light of the cost, and 
effect on the final design. The constraints imposed on the design included a constant blade 
tip diameter dictated by the existing compressor shroud, a pre-determined blade chord set 
by restrictions in the flow path and the need to provide clearance for probe measurements. 
A further requirement was that the stator discharge angle be axial in order to 
accommodate the unique torque measuring system (see Chapter III, Section C2). An 


additional constraint of significant importance was the limitation that the turbine drive unit 


was capable of providing a maximum of 450 to 470 horsepower. 


Within the confines of the physical constraints imposed by the existing compressor 
rig, the designer conducted a preliminary analysis of desired design parameters with the 


goal of producing a stage that achieved, 


“ as high a loading and specific weight flow as was practical, while keeping rotor tip 
Mach number at a moderate level. “ 


Sanger ‘94 
Using the specific head rise parameter, 
WY =C, To (PR -1)/UL 


a parametric study was performed. The specific head rise was calculated over a range of 
tip speeds, and pressure ratios, using as a reference several previous designs, including the 
original Vavra design and the well known NASA Rotor 67. To produce a design that was 
state-of-the-art, the initial selection of design parameters was fairly optimistic. Subsequent 
calculations of the horsepower and inlet ramp angle required to support this selection of 
design parameters revealed that they were in fact too optimistic. The power required 
exceeded that available by nearly 200 hp, and the ramp angle on the spinner hub was 
calculated to be 42 Deg. 

In order to reduce the power required by the stage to that constrained by the 
existing drive turbine, and lower the inlet ramp angle to a more reasonable figure, several 
compromises had to be made. Notable compromises were the reductions in specific weight 
flow and rotor and stage pressure ratio, to alleviate the power imbalance, and a reduction 
in the axial velocity ratio, to reduce the inlet-to-exit area ratio and therefore the inlet ramp 
angle. With these corrections, the required input power was reduced to 457 horsepower, 
and the ramp angle became 28.2 Deg. Selection of additional design parameters such as 
solidity and diffusion factor was primarily based on the designer’s experience with similar 


machines. The results of the parametric study are shown in Figure 3. 
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Figure 3. Results of parametric study [from Ref. 8]. 


3. Final Selection of Preliminary Design Parameters 


Prior to this stage in the design the analysis was restricted to one-dimensional 
methods. To arrive at a detailed description of the flow path, and evaluate the 
performance of the preliminary design, the designer used the AXIDES code. The 
reasoning behind using this code was to reveal any obvious deficiencies in the preliminary 
design and to produce the flow field from which the blading would be derived. The 
AXIDES code provided the full 2-D radial equilibrium solution, where steady, 
axisymmetric flow is assumed, and produced distributions of the flow between blade rows 
at the specified design condition. At this point in the design process the flow-field 
estimation option of the code was used by setting the parameter OP, of the input data set, 


to APPROX (see Appendix A). Using this option, no blades are designed or specified. The 


code simply estimated blade-edge locations from the stacking line, and produces velocity 
diagram and performance information based on these approximations. 

The information that was specified in the input data set to the AXIDES code was 
RPM, mass flow rate, stage pressure ratio, inlet total pressure, inlet total temperature, 
tangential velocity, the flow path of the existing compressor rig, annular station locations, 
and estimated blockage. Although typically input in the form of loss-correlation tables, the 
losses across the blading were implicitly defined by entering the total pressure and total 
temperature distributions at the outlet of the rotor. A discussion of this option of the code 
can be found in Reference 9, page 12. The total-pressure distribution was essentially 
constant except for losses at the hub and tip. The total-temperature distribution was also 
nearly constant, however, higher temperatures were input near the tip to account for 
increased shock losses along the blades. The values used for estimated blockage and total- 
pressure and total-temperature distributions at the outlet of the rotor, were based on 
experience and on the values of successful stages in the same design range. Additional 
information specified in the input data set were the radial distributions of blade leading and 
trailing edge radii, chord length, and maximum thickness-to-chord ratio. All of these 
distributions were based on the designer’s experience with transonic rotors. The input 
used to define the stator attributes were similar. A notable exception was the exit 
tangential velocity which was zero. 

Several iterations through the AXIDES code were conducted, until output 
distributions of such parameters as diffusion factor and loss appeared reasonable. 
Modifications to the input data file consisted mainly of variations in pressure ratio and 
flow path geometry. There was initial concern about the low values of output loss 
distribution across the rotor. It was decided, however, that if an error in efficiency was to 
be made, to design for high efficiency and be wrong was better than designing for lower 
efficiency and be correct. 


At the time of this report only a representative example of the final input data set 


of preliminary design parameters was available. Described as “one of the last” data sets 


used during iterations through the AXIDES code, it was a very close approximation of the 
actual input. This input file was run on the version of the code installed at NPS. The 
output received was the designer’s final preliminary design flow-field and the associated 
performance parameters. Reference 10 contains this input file in its raw form, and the 
associated output file. Appendix A gives a description of input and output file structure 
and parameters. 

Once satisfied with the 2-D flow field and performance, the preliminary design 
parameters were fixed. A list of first cut, and final selection of design parameters is shown 


in Table 1. 


ee 
[TipSpeed SSCS 
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Table 1. First cut, and final selection of design parameters. 


The problem now consisted of blading design to produce the loss distribution, and flow 


angles, specified in the output of the AXIDES code. 


4. Blade Shape Definition 


The AXIDES code provided the designer with a picture of the desired flow field, 
between blade rows, on the meridional plane. The next step was to fit blade shapes to the 
velocity diagrams to produce this flow. From the conception of the design, it had been 


decided that CFD would be used to arrive at these blade shapes. The idea was to make an 


initial “guess” of blade shape, using output from the AXIDES code, run this blade through 


the CFD analysis, and check the results. If the output seemed reasonable the design was 
complete. If not, blade shape modification was required. The designer created initial blade 
sections essentially by hand, guided by the flow angle distributions produced. by the 
AXIDES code and experience with supersonic flow across blading. Iteration on blade 
shape, which turned out to be required, was performed using a NASA in-house blade 
element code, derived from the geometry portion of the AXIDES code. 

The procedure used to screen the initial blade sections was modified over the 
course of the design. Initially, a quasi-three-dimensional code, written by John Denton 
(Ref. 11], was used. This code was coupled with an integral boundary layer code [Ref. 
12], to check the surface boundary layer condition. The Denton code used was an Euler 
code containing a simple transpiration model to estimate the boundary layer blockage. 
Together these codes provided the designer with a quick tool for screening the blade 
shapes. 

During the process of designing the blading, the full three-dimensional version of 
the Denton code (TIP3D), as reported in Denton 1986 [Ref. 13], became available. This 
code was simply an upgrade in the original Euler code, containing a simple approximation 
for viscous effects and an empirical equation defining the distribution of shear stress from 
the wall. Consideration was given to using a code with a more complex and possibly more 
accurate turbulence-modeling scheme, however, given the iterative nature of the design 
process, the increase in accuracy was deemed insignificant when compared with the 
increase in computational execution time. 

Using the quasi-3D results as a “base design”, the iteration on blade shape was 
continued in the same manner as previously described, using the updated version of the 
Denton code. The criteria used for acceptance of the blade included the requirements to 
minimize shock strength in the tip region and to prevent or delay boundary layer 
separation. The shock strength was reduced by minimizing the supersonic acceleration of 
the flow approaching the shock. This was accomplished by making the suction surface at 


the leading-edge portion of the blade wedge shaped, and curving it after the shock to 


10 


provide the required turning. To prevent Separation, the attempt was made to control the 
diffusion towards the trailing edge of the blade. Final-design blade sections, used to 
generate the now-rotating shape for fabrication, labeled “Sanger Design”, are shown in 
Figure 7, located in section B3 of this Chapter. Because these sections were used for 
fabrication, they represented elements of the blade on cylindrical planes, parallel to the axis 


of rotation of the rotor. 


5. Fabrication 


After all criteria had been met, to within certain acceptable tolerances, the 
aerodynamic portion of the design was complete. What remained was the mechanical 
analysis and fabrication of the blading. Historically, the Engineering Design Department at 
NASA-Lewis had used the same geometry specifications as prescribed in the AXIDES 
code, to perform their mechanical analysis. Therefore, the final blade parameters were 
input back into the streamline curvature code using the output option (OP), COORD (see 
Appendix A). This option created an output file containing blade coordinates for 
fabrication. A copy of the input file used to produce the final blade fabrication coordinates 


and the output file containing these coordinates can be found in Reference 10. 


B. DESIGN COMPARISON 


1. Overview 


Streamline curvature methods have been used for the design of many successful 
machines. Historically, the flow-field between blade rows is calculated using these 
methods, a blade shape is arrived at, and the flow within the blade row, where any 
problems are likely to occur, is checked by a more detailed analysis code. The AXIDES 
code was written with the objective of integrating the design and analysis portions of the 
overall design process in a more efficient manner. It produces through-flow output that 


can be directly input into the blade-to-blade analysis codes T-SONIC [Ref. 14], and 


MERIDL [Ref. 15]. The code also provides input features which allow for corrections in 


the blade shape, determined as necessary from the output of these analysis codes, to be 
easily made. The result is a composite through-flow code, capable of both aerodynamic 
and blade design, that can be interfaced with two codes that analyze the flow within the 
blade rows. Currently, the version of the AXIDES code at NPS, is not capable of 
interfacing with MERIDL. An update to the code should include this feature. 

For the reasons mentioned previously, CFD techniques were used throughout for 
the design and analysis of the new Sanger stage. The AXIDES code was not used as a 
blade design tool. It was simply used to describe the flow-field which satisfied the 
requirements dictated by the preliminary design parameters. No blade shapes were input, 
nor specified by the code, and only velocity diagram and performance information was 
output. The initial “guess” of the blade shapes was generated by hand. Although iteration 
on blade definition was performed using a geometry routine extracted from the AXIDES 
code, the blade sections were stacked based on information produced by the CFD analysis, 
and not within the aerodynamic iterations of the streamline curvature method. This section 
contains a comparison between the final Sanger design and a design produced strictly with 
the use of the AXIDES code; that is, the blading specification and aerodynamic 
performance output by AXIDES, when used in a design capacity. The motivation for this 
comparison was two-fold. First, the Crouse 2-D code was to be used for instruction in 
compressor design at NPS. A comparison with the geometry arrived at using advanced 
analysis methods, was a useful exercise of the code. Second, no similar comparison was 
found between a design resulting from relatively simple, and computationally-efficient 
streamline curvature methods, and the geometry obtained using more detailed flow | 


analysis routines, requiring significant investments in computational time. 


2. Method Of Comparison 
The performance of the Sanger design was taken from Sanger, 1994 [Ref. 8]. The 
blade definition parameters and performance distributions, as described previously, were 


set using the Denton TIP3D code. The AXIDES code was used strictly to provide the 


initial flow-field estimate, and produced coordinates of the final blade design in a form that 
could be conveniently used in the mechanical analysis and the subsequent fabrication 
process. During the design review, when the picture of the procedure used in the design 
became clear, there was initial skepticism as to how the final blade design parameters 
could be input into the AXIDES code, simply to provide the desired form of the blade 
geometry in the output. The concern was that the program would modify the blade shape 
as it attempted to converge on a 2-D radial-equilibrium solution. A solution which would 
obviously be, at the least, slightly different than the results of full 3-D CFD calculations. 
This concern was reconciled by checking the blade geometry parameters in the final output 
of the AXIDES code, that was sent to the NASA-Lewis Engineering Design Department, 
against the results of the CFD design, as reported in Sanger ‘94 [Ref. 8]. Specifically, such 
blade parameters as, radial distribution of maximum thickness, and maximum thickness 
location, extracted from the AXIDES output, were reviewed for discrepancies; none were 
found. The blade shape produced by the AXIDES code was identical to the blade 
specified in the input. This appeared to have been accomplished by completely specifying 
all blade parameters in the input file. The aerodynamic Output was of no interest to the 
designer. Aerodynamic input data were used simply to fill in the spaces of the input data 
set and did not represent the performance generated when the code was run. Again, the 
AXIDES code was used simply to reproduce the geometry of a blade designed using CFD 
in the form required by the NASA mechanical design system. 

The design produced using the AXIDES code, and henceforth referred to as the 
“conventional design”, was to represent what the Sanger design would have looked like if 
it had been sent for mechanical analysis directly after the initial design parameters were 
fixed. The blades were designed by taking the final form of these initial parameters, which 
Sanger used to build his “guess” of the blade shape, and inputting them back into the 
AXIDES code. In this case, however, the input file was modified to produce blade shapes 


for fabrication instead of merely calling for flow-field description, as had been done in the 


Sanger design. This was accomplished by changing the OP parameter to COORD and 


using defaults for incidence angle and deviation angle by choosing 2-D for input 
parameters AA, and BB. The 2-D option resulted in the use of calculations specified in 
NASA SP-36 [Ref. 16], for these angles. Additional modifications to the input file 
included setting option CC to OPTIMUM, which assigned a value to the turning rate ratio 
of the blade element segments using an empirical function of inlet relative Mach number; 
setting option DD to shock, which placed the transition point at the location of the shock; 
and setting option EE to TRAN, which located the maximum thickness point at the 
transition point. A copy of the input file, and output file are located in Reference 10. A 
detailed description of the AXIDES code available input options, can be found in 


Appendix A, Section 2, and Reference 9. 


3. Results Of Comparison 


As previously noted, the intent of the Sanger design was to produce blading with 
increased loading at lower tip speeds. A plot of rotor diffusion factor, indicative of loading 
on the blade, is shown in Figure 4. It is clear from this plot that the loading on the Sanger 
rotor is larger over 90% of the blade than that of the “conventional” design. Using the 
detailed picture of the flow-field provided by the CFD code, the designer was able check 
for flow separation along the suction surface of the blade as he increased the amount of 
turning imparted to the flow. Figure 5 shows the amount of energy addition provided by 
the rotor in the form of total temperature ratio. Figure 6 illustrates the increase in total 
pressure ratio across the rotor. The plots indicate an increase in total pressure ratio and 
energy addition produced by careful analysis of the flow within the blade row. The blade 
sections for the conventional design and the Sanger design are shown in Figure 7. This 
figure displays the blade sections at 90%, 50%, and 20% span, from the hub. Looking 
closely at these sections, it can be seen that the curvature on the suction surface of the 
Sanger rotor, near the leading edge, is very small. This was done to minimize the 
supersonic acceleration of the flow as it approached the shock, and therefore reduce its 


strength. The AXIDES code used in the design of the conventional blade also attempted 
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to account for the effects of a large inlet relative Mach number. This was accomplished by 
activating the blade definition option OPTIMUM, located in the input data set. It can be 
seen, however, that the curvature in the leading-edge portion of the conventional blade, 
although modified, is still larger than on the Sanger blade. The turning of the flow was 
accomplished further aft on the blade, as can be seen in the distribution of blade maximum 
thickness location shown in Figure 8. The losses across the rotor were kept nearly the 
same as the conventional blade. A plot of the distribution of loss coefficient is shown in 


Figure 9. 
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Figure 4. Rotor diffusion factor. 
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Figure 5. Total temperature ratio across the rotor. 
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Figure 6. Total pressure ratio across the rotor. 
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Figure 7. Tip, mean-line and hub blade sections (from top). 
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Figure 8. Rotor maximum thickness point location. 
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Figure 9. Rotor loss coefficient distribution. 


The preceding comparisons show that the Sanger blading inputs a larger amount of 
work than blading designed using 2-D methods alone, while maintaining very similar 
through-flow and loss profiles. Assuming that the CFD code used in the design predicts 
the 3-D flow field accurately, the design represents a step forward toward optimizing the 


design of highly loaded, efficient, transonic blading. 
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IH. PREPARATIONS FOR TESTING 


A. OVERVIEW 


The transonic stage designed at NASA-Lewis, was designed to be installed in the 
NPS/TPL transonic compressor test rig, which is shown schematically in Figure 1. Since 
the rig had not been operated since 1986, an evaluation of its condition, and overhaul of 
Some its components was necessary. This section gives an account of the preliminary 
evaluation and subsequent overhaul, to date, required to prepare the machine to test the 


new design. 


B. PRELIMINARY TEST RUNS 


Preliminary test runs of the transonic compressor test rig were conducted in early 
1995. The tests were run with the existing machine as it was left, with the initial rotor, 
designed Dr. Vavra, still installed. This testing was conducted to become familiar with the 
operation of the test rig and to determine the extent of overhaul necessary to return it to 
routine working order. The operating procedure is given in Appendix C. 

A total of four tests were run over a period of two months. The first test was 
considered a basic operational check to make sure the machine would turn over with no 
major component failure. The rig was accelerated up to 5,000 RPM, un-throttled, with the 
inlet piping removed to prevent any deposits in the piping from being pulled into the 


blading. During this test the following mechanical aspects of the machine were monitored, 


1. Pressure ratio across the laboratory 12 stage axial compressor supplying air to 
the drive turbine. 

2. Balance air used to provide relief of the axial load on the high speed 
compressor bearings. 

3. Test rig rotational speed. 


4. Drive turbine and compressor bearing temperatures. 
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All of these indications were monitored at the test control panel. A photograph of the 
control panel is given in Appendix C. 

The results of the operational check were encouraging. Although, the test rig had 
not been operated for many years, there was no appearance of malfunction or unusual 
operation. The tests that followed were conducted in the same manner, however the RPM 
was increased with each run, to a final value of 20,000 RPM. During the testing, the 
bearing temperatures were monitored and compared with values recorded from test runs 
conducted years earlier. These values remained nearly identical to those previously 
recorded until approximately 15,000 RPM (50% design speed). Above 15,000 RPM the 
bearing temperatures increased at a greater rate than previously documented. By adjusting 
the oil-mist cooling air and oil supply settings, the temperatures of the bearings supporting 
the turbine-drive shaft and the downstream end of the compressor shaft were controlled to 
within 2 degrees of earlier recorded values. This was accomplished by reducing the oil 
mist drop rate from approximately 80 drops per minute to 12 drops per minute, and 
increasing air supply pressure from 20 psi to 35 psi. The bearings supporting the upstream 
end of the compressor drive shaft, however, failed to respond to the changes in the 
cooling system. They remained approximately 25 degrees high and continued to increase 
at a rate of approximately 1 degree per hour. 

The decision was made to tear down the compressor in an effort to resolve the 
problem with cooling of the compressor inlet bearings, and to check the bearings for any 
damage that might have occurred due to the high temperatures. Additionally, a decision 
was made to leave the turbine-drive unit intact since no problems in its operation had been 


revealed. 


C. OVERHAUL 


The transonic compressor test rig was a unique machine that had undergone 


several modifications over the years. Although some modifications were documented, the 


overhaul that was conducted for the present project revealed that some of the changes 
made to had not been recorded. Additionally, the procedure involved in disassembling and 
reassembling the test rig was undocumented. The overhaul presented a significant 
challenge because of the mechanical complexity inherent in the arrangement of a high- 
speed air turbine driving a high speed research compressor. The goal was first to 
determine the method by which the rig was to be disassembled and then to carefully 
disassemble and reassemble without causing damage to any of the mostly non-standard 
components. During the course of the tear down and rebuild effort, video tape was taken 


to document the procedure. 


1. Cooling Of The High Speed Bearings 


The driving consideration behind the overhaul was the problem encountered with 
bearing cooling. Disassembly of the compressor hub revealed that, of the two cooling lines 
intended to provide oil-mist air to the compressor inlet bearings, one was completely 
disconnected. A second line which had not been indicated on the original design drawings, 
was bent so that no air could flow through it. Figure 10 shows a drawing of the 
compressor hub assembly. The bearings, notably the ones to the right, were not receiving 
the oil-mist air supply needed to cool them. To correct the problem the cooling lines were 
fitted with new tubing. As an indication of previous problems with bearing temperatures, a 
one-eighth inch copper water cooling line was found encircling the inner compressor hub. 
This water cooling device was not shown on the original drawings, nor recorded in any 


available test log. 


2. Torque Measurement system 


The torque acting on the test rig rotor drive shaft was measured by the deflection 
of a cantilever beam. The beam extended between the outer portion of the compressor 
hub, that rotated on bearings, and the inner hub annulus which was Stationary (see Figure 


10). The rotation of the outer hub was a result of the swirling air-flow exiting the rotor 
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Figure 10. Compressor hub assembly. 
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and impinging on the stator. The stator removed most of the swirl component of the 
velocity, and a honeycomb flow straightener, removed the residual. The stator and flow 
Straightener were both mounted on the outer hub which rotated on the outer hub bearings. 
The cantilever beam (torque balance), whose fixed end was attached to the outer hub and 
free end was constrained by a channel fixed to the inner stationary hub annulus, deflected 
as it tried to prevent relative rotation. Four strain gauges, mounted on the cantilever beam, 
provided the signal that was a measure of the torque imparted to the air-flow by the rotor. 

A static calibration of the torque balance was conducted in preparation for the - 
tests to follow overhaul. A significant drift in the strain gauge output was found. A plot of 
the calibration is shown in Figure 11. The calibration was performed using weights hung 
from a 20.01 inch moment arm which was mounted on the outer hub. The load was 
increased from zero to 35 pounds in five pound increments. A zero of -1.0 milli-volts was 
set on the output voltage display. This was done because the beam appeared to have some 
pre-load on it and the bridge output could not be set to zero. The plot clearly indicates 
repeatability, however, the measurements, conducted for two static load and unload 
cycles, were taken after the Output settled down after a drift of approximately two milli- 
volts per second. This required waiting an average of 15 minutes from the time the load 
was changed until the time the readings were taken. The time dependent nature of the 
Static load, as well as other difficulties with the torque measurement, were documented for 
an early compressor test program using the same rig, by R.P. Shreeve [Ref. 17]. 

After careful consideration of the problem, and the ramifications of possible 
changes to the system, modification of the torque balance assembly was made as shown in | 
Figure 12. The first change was an attempt to insure the fixed end of the cantilever beam 
was firmly “encased” in the outer hub. Four set Screws, mounted in the outer ring were 
installed to apply pressure to the top of the beam and secure it in place. They replaced the 
single (center) set screw in the original design. Additionally, the channel fixed to the inner 


hub annulus was modified. The side of the channel resisting the load was drilled and fitted 
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with a three-eighths inch ball bearing. This was done to insure a point load would be 


transmitted at the free end of the cantilever beam. 
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Figure 11. Torque balance calibration (first). 
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Figure 12. Torque balance modifications. 
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A second calibration of the torque balance was conducted in the same manner as 
the first with the exception that the zero-load output was set at 0.1 milli-volts. A plot of 
the calibration is shown in Figure 13. The results indicated that the drift in the output 
reading had disappeared. However, as can be seen clearly in Figure 13, there was a distinct 
lack of repeatability over the four load and unload cycles. The output seemed to settle out 
after the second cycle, however, the overall uncertainty for the four data sets was in the 
vicinity of 10%. Since the uncertainty in the torque measurement gives an equal 
uncertainty in the measurement of efficiency the lack of repeatability is unexceptable. 
However, since the intent of the overhaul effort was to return the rig to stable mechanical 
operation, no further modifications were attempted. Recommendations for additional 


modifications to the torque assembly include, 


1. Installing additional set screws to attach the outer ring firmly to the outer hub, 
to prevent movement of the parts constraining the fixed end of the cantilever 
beam. 

2. Fitting the existing beam with new strain gauges. 

3. Replacing the existing beam and strain-gauges with a new design having a 


larger, rectangular fixed end. 
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Figure 13. Torque balance calibration (second). 


D. FINAL TEST RUN 


A final test run of the transonic compressor was conducted up to a speed of 
20,000 RPM. This testing was performed un-throttled, with the inlet piping removed. The 
test rig setup was similar to the preliminary test runs except for the addition of a 
honeycomb flow straightener and the original stator that were installed during overhaul. 
The duration of the test was two hours and twenty minutes. The objective of the test run 
was to check the mechanical operation of the rig, and in particular, the behavior of the 
bearing temperatures, following overhaul. 

The test rig operated smoothly during the entire run, at no time indicating 
problems with balancing or vibration. At a speed of 15,000 RPM, the temperature of the 
bearings supporting the inlet end of the compressor drive shaft were stable, but reading 
approximately ten degrees high, relative to temperatures recorded during earlier test runs 
at the same RPM. As the rotational speed was increased from 15,000 RPM, to 20,000 
RPM, the temperature of these bearings increased significantly. The temperature was 
recorded to be 143 Deg. F, and climbing at a rate of approximately five degrees per 
minute, when the RPM indicated 18,000. All other bearing temperatures were stable and 
within limits. The temperature limitations of these bearings require operating levels of less 
than 160 Deg. F. The RPM was reduced to 10,000 RPM and after temperatures stabilized 


the machine was shutdown. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A case study was initiated of the design of a new transonic compressor stage 
which is to be tested at the Turbopropulsion Laboratory at the Naval Postgraduate 
School. The stage was designed by Nelson Sanger of the NASA Lewis Research Center, 
using 3-D, CFD methods. The program of testing, scheduled to follow installation of the 
new stage, will provide data with which to validate those methods used in the design, and 
Serve as a test case for still-emerging analysis codes. 

The process followed by the designer was reviewed and a version of the code 
(AXIDES), which was used early in the design was installed and made operational at the 
School. The code was used successfully both to examine the design process for the new 
Stage, and later to carry out the design of a multi-stage core compressor. In preparation 
for testing, the test rig was operated, partially overhauled, and operated again to 75% 


design speed. The following were concluded: 


Stage Design 


e The initial flow-field used by the designer was reproduced successfully using AXIDES, 
and the (conventional) blading associated with this flow-field was generated. 

¢ More highly loaded blading was produced using CFD advanced design methods than 
could have resulted from the use of the AXIDES code alone. 

¢ The differences in blade shape that enabled the higher loading were, visually, relatively 
subtle. 

e The new stage was manufactured successfully and is ready to test after installation in 


the nig. 


3] 


AXIDES Code 


® The AXIDES code was made operational on NPS, AA Department workstations, and 
on the TPL PC. 
® The code was run for several design cases, the input and output of which serve as 


examples of the varied uses of the code. 


Test Ri 


e Preliminary test runs were conducted up to 20,000 RPM (design speed is 30,000 
RPM). 

© The machine was overhauled to correct problems with cooling of the compressor inlet 
bearings. The bearings were in good condition and were re-installed. The cooling air 
supply system was checked and returned to working order. 

e The torque measuring system was modified; however, a lack of repeatability remains. 

e Following the reassembly, a test run was completed up to 20,000 RPM, and indicated 
continuing problems with the bearings at the inlet end of the compressor. 

e The pressure at the rotor exit is greater when the machine is run without the inlet 
piping, than when the piping, screens and throttle are installed. Since the rotor exit 1S 
where the cooling air is exhausted after passing through the bearings, the higher 


pressure may be preventing the cooling air from circulating as intended. 


32 


The following recommendations are made: 


The “conventional” blading design should be put through CFD analysis, using a more 
recent 3-D code, and compared with the results of the Denton, TIP3D code. 

The tools for complete axial compressor design, including preliminary 1-D, 2-D, and 
3-D CFD, codes, all now available at NPS, should be integrated. 

The AXIDES code should be upgraded to include a graphics package, and reliable 
interfacing with the off-design code AXIOFF. 

The transonic test-rig inlet pipe and housing should be installed, and the test rig 
operated, to determine if the reduced pressures at the rotor exit effect the compressor 
inlet bearing temperatures. oe 

The new stage testing program should move ahead without delay to prevent the recent 


experience from being lost. 
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APPENDIX A. AXIDES CODE 


A. OVERVIEW 


The AXIDES code written by Jim Crouse is a composite aerodynamic and blade 
design code. It uses the streamline curvature method to arrive at a full, 2-D, radial 
equilibrium solution of the flow between blade rows in the meridional plane, assuming 
steady, axisymmetric flow. The code can be used for either aerodynamic or blade design, 
or both. In the case where both design capabilities are required, the code performs four 
iterations of the aerodynamic solution, and then stacks the blades. This process is 
continued until the solution converges. The final blade shape is calculated after the 
aerodynamic solution is complete. The code can be used for either single or multistage 
axial compressors. A complete description of the AXIDES code can be found in 
Reference 9. 

A version of the AXIDES code was purchased by the Naval Postgraduate School. 
The code is originally PC based and has been modified to run on the Department of 
Aeronautics and Astronautics workstations. The PC version was installed and aneried on 
the Turbopropulsion laboratory’s 486 machine. After Overcoming initial difficulties with 
input and output assignments and initialization of variables, the code was used for the 


following, 


1. To reproduce the initial flow-field solution used in the design of the new, NPS 
transonic compressor stage. 

2. To reproduce the geometry of the final blade design, of the same stage, used in 
the mechanical analysis and fabrication of the hardware. 

3. To design a stage using the strictly 2-D calculations of the code, for the 
comparison with the new stage, arrived at using 3-D CFD methods. 

4. To design a multi-stage compressor for the AA Department Engine Design 


Course. 
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B. INPUT 


Creating the input file to the AXIDES code is crucial to its successful use as a 
compressor design tool. Along with the need for a basic understanding of internal and 
turbo-machinery flows, a general understanding of the input file structure is required to 
produce realistic output. This section provides the tools necessary to build a working input 
file. 

The input to the code is in the form of a formatted data field. Each piece of the 
data has a specific location in the data field and a specific length. If the input data is not in 
the correct location, the formatted READ statements in the code will not read the correct 


information and the program will run based on invalid input, or may not run at all. 


Note: 


An important note here is that a “dummy” editor, such as Notepad’™ on the PC, or 
jot, on the workstations should be used to create or modify the input file. Other 
editors including DOS and word processing editors have been found to change the 
formatted structure of the file, however slightly, enough to move the input data out 


of position and prevent the correct input from being read. 


Figure 14, located on page 49, shows an example working input data file. This file 
was provided with the code in an effort to get the program up and running on the AA 
Department machines. To arrive at a successful input file it is recommended that an 
example file be run first, to become familiar with the operation of the code, and then 
modifications to that file be made to create an original design. Five input files were used to 


produce this report. Four of these, namely, 
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e Examp.dat : The example file provided when the code was purchased, 

e Approx.dat : The file used to set final preliminary design parameters, and 
define the desired flow-field for the new Sanger stage, 

e Conv.dat : The file used to produce the “conventional” stage for the 
design comparison, 


¢ Multax.dat : The file used in the design of a four stage compressor, 


are provided for use as a starting point to build an input data set. The input file Examp.dat 
is shown in Figure 14. A copy of the output associated with Examp.dat (Examp.out), is 
located in Appendix B. The last three files listed above can be found in Reference 10. The 
fifth input file, also located in Reference 10, titled Sanger.dat, was used simply to produce 
fabrication coordinates and the aerodynamic input is not representative of the blades that 
the code designed; See Chapter II, Section B2, paragraph 1, for a discussion of this. 

The input file to the code is broken up into two main sections. The first section 
contains general information which includes the number of stages, the flow path or casing 
and hub geometric coordinates, the RPM, and the desired overall pressure ratio. The 
second section consists of annular calculation station and blade row data. This data 
includes axial locations of annular stations and blade rows, boundary layer blockage 
estimates, the number of blades on the rotor and Stator, solidity, and blade definition 
parameters such as leading and trailing edge radii. 

Figures 15-19, located on pages 50-54, provide another look at the example input 
data file mentioned previously (Examp.dat). This illustration is included as an aid in 
identifying the input file parameters used for this particular design. It must be emphasized 
that this figure is not a working file but is included only for parameter clarification. The 
corresponding working input file is shown in Figure 14. Also included is an illustration of 
the general form of the input (Figures 20, 21, and 22, pages 55, 56, and 57), taken from 
notes provided by the author of the code, with all possible input parameters and associated 


formatting. For a complete description of the code and its input see Reference 9. 
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The following is a description of the input parameters used in the example input 


file Examp.dat. It addresses each parameter in the file in the order shown in Figures 13-17. 


The information provided below was extracted entirely from Reference 9, and is included 


merely for the readers convenience. 


e | 

e IROTOR 
e J 

e K 

e Column 
e TITLE 

e ITG 


Indices 


Calculation station index 
Rotor index. 
Streamline index. Numbered from one at the tip. 


Loss set index. 


General Input 


Not an input parameter. The first three lines simply show the 
relationship of the data below to the 80 column FORTRAN data 
field. 

The title of the design. 

Output blade fabrication coordinate specification. The output of the 
blade shapes produced by the code is in the form of suction and 
pressure surface coordinates on blade sections. These coordinates 
are located at uniform and round number increments of the 
chord-wise distance from the leading-edge. Included in the table of 
coordinates are the tangency points between the blade surfaces and 
the leading and trailing-edge ellipses. The locations of the point of 
tangency for the suction and pressure surfaces will be different. If 
the parameter ITG is set to 0, the program will insert an 


artificially high value of 99.9999, for the height of the surface 
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ITG (cont’d) 


IBR 


NSTRM 
NROW 


opposite the surface that is tangent at that chord-wise distance. This 
Clearly identifies the junction point between that surfaces and the 
end ellipse. The leading and trailing-edge information is provided 
below the table of coordinates. If ITG is set to 1, the program 
simply outputs the actual height of the surface opposite the surface 
that is tangent at that point. The tangency points must then be 
determined by noting the change in the surface heights near the 
leading and trailing-edges. 

Set equal to 0 to get error messages printed when running the 
code. Set equal to 1 to suppress the error message output. If set 

to a value greater than 1, the program will provide no error 
messages during the run, and stop after the performance summary 
is printed, printing neither radial distribution information, nor 
fabrication coordinates, even if specified in the input file. 

Number of streamlines (11 max.). 

Number of blade rows (20 max.). 

Number of annular stations. Fifty maximum, including leading and 
trailing edges of each blade. Must be at least four prior to first blade 
row and three after the last. 

Number of loss sets input in the form of tables of correlation data. 
Number of outer case wall geometric coordinates. 

Number of hub geometric coordinates. 

Compressor RPM. 

Mass flow rate (Ib./sec). 

Desired overall pressure ratio. 


Molecular weight of air. 
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SCALEF 


CP(K) 
LUNITI 


LUNITO 
FLO(K) 
TO(1,J) 
PO(1,J) 
VTH(1,J) 
XTIP(I) 
RTIP(I) 
XHUB(I) 
RHUB(I) 


DLOSS(K,J,1) 
DFTAB(K,J,1) 


Scaling factor. If set to other than 1.0, the program will scale all 
linear dimensions and the weight-flow to a different size 
compressor, by the factor indicated. 

Constants for specific heat polynomial. 

Units of input. Set equal to 1 if SI and 2 if English. 
Desired units of output. Set to 1 for SI and 2 for English. 
Cumulative flow fraction at streamlines, from the tip. 
Inlet total temperature at streamlines from tip. 

Inlet total pressure at streamlines from tip. 

Inlet tangential velocity at streamlines from tip. 

Axial coordinate of casing wall points. 

Radial coordinates of casing wall points. 

Axial coordinates of hub wall points. 

Radial coordinates of hub wall points. — 

Array of loss parameter correlations (five sets max.). 


Array of diffusion factor correlations (five sets max.). 
Calculation Station And Blade Row Data 


Type of axial station : 
: ANNULAR 
: ROTOR 
- STATOR 
,or, depending on location in data set, 


Incidence angle input 


- 2-D NASA SP-36 correlations performed internal to 
code. 
: 3-D NASA SP-36 correlations. 
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AA (cont'd) 


ZTIP(I) 
ZHUB(1) 
BT (1) 


BH(I) 
BLEED(I) 
DLIM(IROW) 


ALIM(IROW) 


CRENGY 


BMATL(IROW) 
NXCUT(IROW) 


: SUCTION Zero inc. to suction surface. 
: TABLE Input in tabular form as INC(IROW,J), at the end of 
associated blade rows input 
data set. 
Tip axial coordinate of annular station or blade stacking line. 
Hub axial coordinate of annular station or blade stacking line. 
Tip boundary layer blockage for annular station or blade LE or TE. 
Can be input as a fraction of annular area, or as the displacement 
from the wall in inches if preceded by a negative sign. 
Hub boundary layer blockage (same as BT except applied at hub). 
Fraction of weight flow bled off at annular station or blade row. 
Diffusion factor limit. Applies to tip for rotors and hub for stators. 
If limit exceeded the program will reduce the energy addition across 
that particular blade row and try to make it up in others that 
have not exceeded their limit. If all blade rows are at there 
limit the overall pressure ratio will be reduced. 
Minimum relative flow angle limit leaving rotor hub, or, maximum 
allowable Mach number entering the stator hub. Adjustments made 
to satisfy limitations are the same as for DLIM except applied to 
flow angle and Mach number. 
Cumulative fraction of energy input across rotor to that input 
across stage. If greater than 2, interpreted as rotor tip exit total 
temperature in degrees Rankine, and is converted to energy 
addition internal to program using the total temperature profile 
input in the form of parameters PRA-PRE. 
Material density of rotor (Ib./in.’). 
Number of blade sections for which fabrication coordinates will be 


produced. If zero, the program will select this number based on 


4] 


e NXCUT (IROW) aspect ratio. If negative, represents the number of sections desired, 


(cont’d) 


ILOSS 


OP 


BB 


and program reads specific locations in a table input as parameter 

XCUT, at end of entire data file. 

Loss set used for blade row. If less than or equal to 0 total 

pressure at the blade row exit is input instead of losses being 

computed internal to the program. See Reference 9 for a discussion 
of parameters PTT and PTC. 

Input option controlling amount and type of output information. 

- APPROX Only velocity diagram output based on estimated 
blade edge locations. 

- VEL. DIA. Only velocity diagram output based on blade edge 
locations that are input as ZTEMP, and RTEMP, at 
the end of the associated blade rows input data set. 

- DESIGN Output consists of velocity diagram information 
only based on stacked blade edge locations 
computed internal to code. 

- COORD Output includes velocity diagram and blade 
section fabrication coordinate information 
based on stacked blade edge locations. 

Part of incidence angle TABLE input option. If SS, the code will 

reference incidence angle input to the suction surface at the leading 

edge. Otherwise it will be referenced to the leading edge center- 
line. 

Deviation angle input. 

. 2-D NASA SP-36 correlations. 

. 3-D NASA SP-36 correlations. 

: TABLE Input in tabular form as parameter DEV(IROW,,J), 


at the end of associated blade rows input data set. 
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e BB (cont’d) : CARTER _ Deviation calculated internal to the code 


using Carter’s rule. 


e CC Blade element geometry input. 
; CIRCULAR Code Designs circular arc blade sections. 
: OPTIMUM Curvature at the leading edge set by an 


empirical function of inlet relative Mach 
number. Below M’1 of 0.8 the shape will be 
circular arc. As the relative Mach Number 
increases, the ratio of turning rates of the 
front and rear blade segments is reduced to 
avoid large shock losses. 

: TABLE Ratio of front to rear segment turning rates. 
Input in tabular form as PHI(IROW), at the 
end of associated blade rows input data set. 

e DD Input location of transition point. 

: CIRCULAR Transition point between front and rear 
segments put at mid-chord. 

>: SHOCK Transition point located at suction surface 
shock attachment point. 

: TABLE Input as TRAN S(IROW,J), in tabular form, 


at the end of the associated blade rows input 


data set. 
e EE Maximum thickness location input. 
: TRAN Locates the max. thickness point at the 


transition point. 
: TABLE Maximum thickness location input as 
ZMAX(IROW,J), in tabular form, at the end 


of the associated blade rows input data set. 
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e EB 


¢ CHOKE(IROW) 


e BLADES(IROW) 
e SOLID(IROW) 

e TILTUROW) 

e PRA-PRE(IROW) 


¢ TALE-TDLE(IROW) 


e© TATE-TDTE(IROW) 


e TAMAX-TDMAX(IROW) 


Used with TABLE option of max. thickness 
location. If LE, the maximum thickness location is 
input as a fraction of blade chord from the leading 
edge. Otherwise it is input as a fraction of the chord 
from the transition point. 

Choke margin input. If greater than 0 the program 
will adjust incidence angle in an attempt to provide 
this margin. If 0 the code makes no adjustment. 
Number of rotor or stator blades. 

Input solidity. 

Stacking axis tilt angle. 

Coefficients of the input pressure or temperature 
profile polynomial at the exit of a rotor, or 
tangential velocity at the exit of a stator. The 
pressure profile is input here if stage energy addition 
is input using CRENGY. If the rotor tip exit total 
temperature profile is input as the parameter 
CRENGY, the total temperature profile polynomial 
coefficients are input here and the pressure profile 1s 
input using PTT and PTC. 

Polynomial coefficients describing the distribution of 
the ratio of blade element leading-edge radius to 
chord. 

Same as TALE-TDLE except applied to the trailing- 
edge. 

Polynomial coefficients describing the distribution of 


the ratio of blade element maximum thickness to 


chord. 
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e CHORDA-CHORDC(IROW) Polynomial coefficients describing the distribution of 
the ratio of blade element chord to tip chord on a 
projected plane. | 

e IDEF(IROW) Blade section definition parameter. If 0 the blade 
element surfaces and centerline are described by 
dk/ds = const. If not 0 they are defined by a fourth- 
degree polynomial input on the next line of the data 


set. See Reference 9, for a more complete 


discussion. 
e INC(IROW,)) Incidence angle array on streamlines. 
e DEV(ROW,J) Deviation angle array on streamlines. 
e PHI(IROW,J) Inlet/outlet segment turning rates on streamlines. 
e TRANS(ROW,J) Transition point location‘on streamlines. 
e ZMAX(IROW,J) Maximum thickness point location on streamlines. 


C. CODE EXECUTION 


The name of the executable file running on the AA Department workstations, and 
a 486 PC at the Turbopropulsion laboratory is Stream. To run the code on the 
workstations simply type stream and enter. When running the code on the PC, you must 
first activate the extended memory by typing 0s386 and enter, this is a requirement of the 
FORTRAN compiler, and then up stream and enter. At this point you will be prompted 
for the name of the input data file. This data file must be located in the same directory as 
the executable file. If you stick to a prefix of less than seven alpha-numeric positions, and 
a suffix of three letters, e.g. NASAR20.DAT, you will be O.K. After typing in the name of 
the input file press enter and you will immediately be prompted for the name of the output 
file. Use the same direction for the name of this file as was used for the input file. After 
typing in the output file name press enter again. If the code has any trouble reading the 


input file the program will stop at the point where the difficulty was encountered. By 
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checking the output file that, will be printed up until the time of the error, you should be 


able to find the source of the error in the input. 
Reminder : 


An important note here is that a “dummy” editor, such as Notepad’™ ,on the PC, or 
jot, on the workstations should be used to create or modify the input file. Other 
editors including DOS and word processing editors have been found to change the 
formatted structure of the file, however slightly, enough to move the input data out 


of position and prevent the correct input from being read. 


Once the input is read the program will ask if you want an off design output file to 
be created. At the time of this report the off-design code that was written to interface with 
the design code was undergoing modifications. If yes or Y is selected the program will 
create an input data set to the off-design code, however, given that the off-design program 
is not up and running it would be prudent to type N. After the off-design prompt the 
program will enter its aerodynamic and blade design (if activated) iteration. It is 
recommended that the input parameter IBR be set to zero until the program is running to 
take advantage of any error messages displayed during these iterations. After the programs 
aerodynamic solution converges, you will be prompted to reset blade angles. This option 
has not been exercised. It is up to the reader to investigate this “uncharted territory”. If N 
is entered at this prompt, the program will continue by stacking the blades one final time | 
and printing the coordinates. The output and, if selected, the off-design input file will be 


placed in the same directory as the executable and the input files. 
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D. OUTPUT 


The output of the code is broken up into four sections, consisting of the following; 


1. Asummary of the input. 

2. Velocity diagrams and flow-field information at annular and blade edge 
calculation stations. 

3. Stage performance summary. 


4. Fabrication coordinates if blade design option is activated. 


Additionally, iteration output is provided, after the input summary, if the error message 
parameter IBR is set to zero. A copy of the output for the five runs of the code, as 
described in part A of this appendix, is located in Reference 10. 

The output of the code is formatted in a landscape orientation. Because of the 
length of the file, up to more than eighty pages for a multistage compressor, the file must 
be printed from a word processor on a PC. The conversion process involved in 
transferring the file from DOS to the word processor tends to change the layout. Attempts 
to print the correct layout of the output on a PC were unsuccessful. However, using the 


following command, 
>/usrfocal/bin/a2ps.5.2 -p -nH -l -1 -ns -nP -F8.5 FILENAME | Ip -dhp3si_1 


on the AA Department workstations, it was possible to print the output. The current » 
operation of the networking system requires this exact input to specify all the parameters 
needed for printing. 

A final note concerning the output is that if the program is run using the PC 
version, the output is broken up into two parts. The first part uses the name defined by the 
user, and contains the performance information. The second part is located in a separate 


file, using the same prefix and .COO for the suffix, and contains the blade fabrication 
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coordinate data. The fabrication coordinates of the output blade design are given for the 
suction and pressure surfaces on planes parallel to the axis of rotation, relative to the 
chord, and relative to the axial direction (“turbomachinery orientation”). When running the 
program on the workstations these files are combined in the user defined output file. 
Recall that the fabrication-coordinate data are only output when the blade design option is 


activated. 


E. CONCLUDING REMARKS 


As noted previously, the AXIDES code is a 2-D code that produces the 
aerodynamic and blade design (if desired), of single or multistage axial compressors. The 
run time of the code is measured in seconds, and the output of the code is in a form that 
can be readily used in mechanical analysis routines and in the fabrication process. The 
AXIDES code is available commercially and can be conveniently used on a personal 
computer. Although the code is based on 2-dimensional flow theory, if the application for 
the compressor being designed does not require “cutting edge” design technology, the 
short run time and conveniently formatted output make the AXIDES code a powerful tool 
to use. It is understood that smaller companies that design compressors for use primarily 
in an industrial capacity, and not for high performance aircraft, currently use the AXIDES 
code as a final design tool. These companies can not make a major investment in the 
development of their own design system and “reverse engineering” data base as do the 
larger aircraft engine companies. Equally, they can not afford the engineering and 


computer costs of using sophisticated CFD analysis codes. 
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Figure 14, Examp.dat input data file. 
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rarer renner ne ee ee ee ee eile COLUMN ee REEERE KEE KERE KERR K EK EEKEEWeE HK 


£1111111112222222222333333333344444444445555555555666666666677777777778 
1345678901234567890123456789012345678901234567890123456 7890123456 7890123456 7890 


TITLE ITG IBR 

First Stage Redesign AF WASH 2 Stage - AR =1.52 a % 

NSTRM NROW NA NLOSS NTIP NHUB ROT FLOW PRATIO MOLWT SCALEF 
oo a 1.0000 
CP(K).,<K=1,3 LUNITI 

.23762070E+00 qeecesson i - .28462996E-06 se 

CP(K) , K=4,6 LUNITO 

.81650840E- 09 Lonessaeee- .28442579E-15 a 


CUMULATIVE FLOW FRACTION AT STREAMLINES 


ne ee 
ow ee Serr S SS SSS SSS SS SS eee eee 


.1257 .2457 .3600 .4686 .5714 . 6686 . 7600 . 8457 
.9257 1.0000 


a 
a — i ened 


288.17 288.17 288.17 288.17 288.17 288.17 288.17 288.17 
288.17 288.17 288.17 


INLET TOTAL PRESSURE AT STREAMLINES 


—_——— a 
I eee 


10130.0 16130.0 10130.90 10130.0 10130.0 10130.0 10130.0 10130.0 
10130.0 10130.0 10130.0 


INLET TANGENTIAL VELOCITY AT STREAMLINES 


nn ee 
weer rere sesstessssesseee=5-5-5>--55°-5e0e eee 


.0000 .9000 .0000 .0000 .0000 .0000 .0000 .0000 
.0000 .0000 .0000 


a 


-35.5600 -30.4800 -25.4000 -20.3200 -15.2400 -12.7000 -10.1600 ~7.6200 
-5.0800 -2.5400 0.0000 4.0640 5.0800 7.6200 10.9730 12.7000 
15.2400 17.7800 20.3200 24.0000 


Figure 15. Examp.dat input data file (breakdown). 


50 


TKK KK NRE EKA eee Kee COLUMN ie hee ee ee ee ee ee rears 


11111211112222222222333333333344444444445555555555666665666677777777778 
12345678591234567890123456789012345678901234567890123456789012345678501234567890 


RADIAL COORDINATE OF CASING WALL POINTS 


SWS STR SSR SSS SS SBS SSS BZ e tte eee eM ew ee ee Le 
ee ee ee a a ee ee ae oe! oe oe oe a a a 


25.4000 25.4000 25.4000 25.4000 25.4000 25.4000 25.4000 25.4000 
25.4000 25.4000 22.3900 25.1460 25.1460 25.1460 25.1460 25.1460 
25.1460 25.1460 25.1460 25.1460 


AXIAL COORDINATE OF HUB WALL POINTS 


5 i an 
—=_—_———— eS See SS See eS Se See EE SSS SSS Ssse 


-35.5600 -30.4800 -25.4000 -20.3200 -15.2400 -12.7000 -10.1600 - 7.6200 
-5.0800 -2.5400 0.0000 4.0640 5.0800 7.6200 10.9730 12.7000 
15.2400 17.7800 20.3200 24.0000 


RADIAL COORDINATE OF HUB WALL POINTS 


pa 
—_—_—_——— ee ES eS Se SE wee ew eS SS LSS eee eae SST SEZ ZTE 


7.3020 7.3020 7.4930 8.2550 9.1440 9.5890 10.0960 10.5150 
11.0490 11.7470 12.6680 13.9700 14.3360 14.7320 15.2400 15.2400 
15.2400 15.2400 15.2400 15.2400 


LOSS PARAMETER ARRAY DIFFUSION FACTOR PARAMETER ARRAY 

SET #1 
.0167 .0199 .0243 0312 .0406 . 3000 . 4000 . 5000 .6000 . 7000 
-0123 .0143 .0176 .0222 .0289 . 3000 . 4000 .5000 -6000 . 7000 
-9100 .0113 .0132 .0163 -0210 . 3000 -4000 .5000 -6000 . 7000 
.9080 .0089 .0103 .0130 .0165 . 3000 .4000 - 5000 - 6000 - 7000 
.0080 .0089 .0103 -0130 -0165 . 3000 . 4000 .5000 - 6000 . 7000 
. 0080 .0089 .0103 .0130 .0165 . 3000 - 4000 .5000 -6000 . 7000 
.0080 .0089 .0103 -0130 .0165 . 3000 .4000 .5000 .6000 . 7000 
-0080 .0089 .0103 0130 .0165 . 3000 .4000 .5000 . 6000 . 7000 
.0090 .0103 .0122 .0153 .0200 . 3000 -4000 -5000 .6000 . 7000 
.0092 -0110 .0140 .0182 .0243 . 3000 . 4000 .5000 .6000 . 7000 
.0104 .0127 -0168 -0221 . 0296 . 3000 - 4000 . 5000 . 6000 . 7000 

SET #2 

=SeS=S== 
0309 .0336 -0373 . 9430 .0508 . 3000 .4000 .5000 . 6000 . 7000 
.0272 .0290 .0320 .0362 . 0423 . 3000 - 4000 -5000 . 6000 . 7000 
-0250 .0263 0282 0313 .0360 - 3000 - 4000 .5000 . 6000 . 7000 
. 0230 .0239 .0253 .0280 .0310 . 3000 -4000 -5000 - 6000 . 7000 
.0211 .0220 -0234 .0261 - 0296 . 3000 -4000 .5000 . 6060 - 7000 
.0212 .0222 . 0236 . 0264 .0299 . 3000 -4000 .5000 - 6000 . 7000 
.0214 .0226 .0241 . 0269 . 0306 - 3000 -4000 .5000 - 6000 . 7000 


-0218 .0231 .0248 .0278 .0317 . 3000 . 4000 .5000 . 6000 - 7000 


Figure 16. Examp.dat input breakdown (cont'd). 
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Porro l COLUMN rarer aren an enn ne al 


LLL111111122222222223333333333444444449445555555555666666656577777777778 
123456789012345678901234567890123456789012345678901234567890123456 78901234567890 


SET #2 (CONT'D) 


—=wremsartertreaee eee eee 


<O233 0248 202-70 .0303 .0347 . 3000 .4000 5000 .€500 , 1000 

S02 72 “0290 .0320 ‘9362 . 0423 . 3000 .4000 .5000 . 6000 . 7000 

.0294 ,O3177 .0358 .0441 .0486 . 3000 .4000 *5006 . 6000 . 7000 
SET #3 

20 “240 

me) .05 

.0 ee 

.0 

.O 

m0) 

0 

.O 

.9 

me) 

.0 


CALCULATION STATION AND BLADE ROW DATA 


——_ we ee ee 
SerrrtotraesstsoeseSe ssa rT See eee eer er eee 


AA ZTIP 2HUB BT BH BLEED 
ANNULAR -32.0000 -32.0000 0000 0000 0000 
ANNULAR -21.0000 -21.0000 0000 0000 0000 
ANNULAR -13.0000 -13.0000 0000 0000 0000 
ANNULAR -7.90000 -7.0000 0000 0020 0000 
ANNULAR -3.0000 -3.0000 0030 0030 0000 
AA ZTIP ZHUB BT (LE) BH(LE) BLEED(LE) 

ROTOR 2.0500 2.0500 0070 0070 0000 

DLIM ALIM BT (TE) BH(TE) BLEED(TE) CRENGY BMATL NXCUT 

-4600 -20.000 .0100 .0100 .0000 1.0000 7.823 0 

ILOSS OP AA AB BB cc DD EE EB CHOKE 

1 COORD TABLE SS TABLE TABLE TABLE TABLE LE .0000 


Figure 17. Examp.dat input breakdown (cont'd). 


a2 


i i ti hh hee ee eee eee ee ee COLUMN ieee ee eee ee ee ee eres 


11131111112222222222333333333344444444445555555555666666666677777777778 
12345678901234567890123456789012345678901234567890123456789012345678901224567890 


# BLADES SOLIDITY TILT PRA PRB PRC PRD PRE 
44.0000 1.3000 © e900,» «g000, Ss sg000S«000 Swot sang 
TALE TBLE TCLE TDLE TATE TBTE TCTE TDTE 
0055 0000 ©—«.0000 «0.0000 -«.g0ss.«=Ss«w000«S«w000Ss« ano 
TAMAX TBMAX TCMAX TDMAX CHORDA CHORDB CHORDC IDEF 
0327 0496 ~—«.0000 «0.0000 += 0000S o000s«=000 


INCIDENCE ANGLE ARRAY ON STREAMLINES 


i re I 
—— oe eS Ee ee ee Tre wee Se SA ZZ 


.00 . 00 .00 . 00 ~.00 . 00 .00 . 00 
.00 .00 . 00 


eS — SE 
—_—_—_———— SS SS SSS eS Ee eee ewe EES SESS 


SOS SESS SPSS SPS STM SSeS Set Sse ee ee ee ee ee eee 
—_———— Sw HR MOSS SSS SIMS SBS SSH HSS BSS cts Sess Sets 


1.000 1.000 1.000 


ee 
—_————— eS SSK SSS SSSR SIS SSSR SSS StH SSS Sse 


0-.700 0.680 0.660 0.640 0.620 0.600 0.540 0.480 
0.420 0.360 0.300 


MAXIMUM THICKNESS POINT LOCATION ON STREAMLINES 


0.500 0.500 0.500 0.500 0.500 0.500 0.500 0.500 
0.5000 0.5000 0.5000 


AA ZTIP ZHUB BT (LE) BH (LE) BLEED (LE) 


— =< — ——— te === oe ee ee oe oe 
=Ssc= =sS=s= = SSS = ==. SS =-sSs ==> S=Ss==> 


STATOR 9.4000 9.4000 .0100 -0100 . 0000 


Figure 18. Examp.dat input breakdown (cont'd). 
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CREEK EEK ENE KEE K KKK Ke hee COLUMN ara ne ee ee ee ee ee lll 


L1L11111112222222222333333333344444444445555555555666666666677777777778 
12345678901234567890123456789012345678901234567890123456789012345678901234567890 


DLIM ALIM BT (TE) BH (TE) BLEED (TE) NXCUT 
yo00«1.0000=SS~—S0200=SS«—=0100« 0000 9 
ILOSS OP AA AB BB CC DD EE EB CHOKE 
Tt goo “Gaeneesaanie: Gasce Guock> «| “tnenee “0000 
# BLADES SOLIDITY TILT PRA PRB PRC PRD PRE 
nueecsy "S"SJoo. «oon ==«w000«S=«00 «0000S «000.0000 
TALE TBLE TCLE TDLE TATE TBTE TCLs TDTE 
see cosa ~«wo00d-=«0000«=S«S«0E~—«0.0000 «0000-0000 
TAMAX TBMAX TCMAX TDMAX CHORDA CHORDB CHORDC IDEF 
Se. eile. des, odes» geeon  cooon, cooae. “6 


a ee 


0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 


a 
Serta toes sense Se ae Se SS Se Se See eee eee ere eo 


ee 
Sore stem em tH SSS STS ST SS STS SS STS Se SS ee Se eee eer eee 


a 
ar NA i eens 


50 .59 .50 .50 .50 . 50 :59 50 
.50 .50 .50 
AA ZTIP ZHUB BT BH BLEED 
ANNULAR 14.0000 14.0000 .0150 .0150 .0000 
ANNULAR 18.0000 18.0000 .0150 .0150 .0000 
ANNULAR 23.0000 23.0000 .0150 .0150 .0000 


Figure 19. Examp.dat input breakdown (cont'd). 
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Input Data Field for the Program AXIDES 


(General or Overall Data) 


[T] | 
Jy |Stetetetesncacanaewensanccenescenses Colum SEER Lea ete een oe ee ween 
Option ip| TUVVTTI1992222222222333333333346666666466555555555 5 566666666677777777778 | Format 


fe | 12345678901234567890 1234567890 12345678901234567890 123456 7890123456 7890123456 7860} 


7G 18R] 1BA4, 214 


\c| —oTETLECK), K=1,18 
JG|NSTRM WROW NA NLOSS NIIP WHUB ROT FLOW PRATIO. —- MOLWI SCALEF | 615, 5F10.4 

[c| (CP(K),K=1,3) LUNITI[ 3€20.8, 10x, 110 
jc] (CPCK) ,K=4,6) LUNITO} 3E20.8, 10x, 110 


(FLOCK), K=2,NSTRM) | BFI0.4 


JG} (Cumulatative Flow Frac. From Tip & Streamlines) 
(FLOCK), K=10,NSTRM)| 8F10.4 


IF NSTRM>D jc] 


fc] CInlet Totat Temp. & Streamt ines) (TOC1, J), J=1,NSTRA) | 8F10.4 

IF WSTRM>B {c{ (TOC 1,5), 229, NSTRM) | 8F10.4 
|G] CIniet Total Press. 3 Streamt ines) (POC1, J), J=1,NSTRM) | 8F10.4 

If NSTRM>B {G{ (POC 1,3), J=9,NSTRM) | 8FI0.4 
IG} Cinlet Tang. Vel. 3 Streamtines) CVTHC1, 3), JE1,NSTRM) | 8F10.4 

IF WSTRM>B |G | (VTH(4, J), 529, NSTRM) | 8F10.4 
|G] CAmiat Coord. of Casing Watt Pts.)  (XTHEC1), 1=1,NTIP) | 8F10.4 

IF WTIP>8 = {G| (XTIPCL), F=9,NTIP) | 8F10.4 
IF HTIP>16 [G| (XTIPCI), 1217, NTIP) | 8F10.4 
IF WTEP>24 [GC] CXTIPCI), 2225, NTIP) | 8F10.4 
IF WTIP>32 [Gf CXTIPC(L) Pe33,NTIP) | 8F10.4 
[G} (Radial Coord. of Casing Wall Pts.) (RTIPCE),1=1,NTIP) {| 8F10.4 

IF NTEP>8 = |G| CRTIPCI), L=9,NTIP) { 8FI10.4 
IF WTEP>16 [Gf CRTIPC1),Pet7,NTEP) | 8F10.4 
IF NIIP>24 {G] CRTIPCT), 1825, NTIP) | 8F10.4 
IF WTIP>32  {G| CRIIPCI), F233, NTIP) | 8F10.4 
|G{  CAxial Coord. of Hub Well Pts.) (XRUB(T), Det, WTIP) { 8FI0.4 

IF wHUB>B = |G] (XHUB(T), F=9, NTIP) f 8F10.4 
If WHUB>16  [G] (XHUB(T), P17, NTIP) | 8F10.4 
If WHUB>24 |G! (XHUB(1), F225 ,NTIP) } 8F10.4 
If WHUB>32 = |G| (XHUB(T), 1 =33,NTIP) | BFI0.4 
{G]  (Redial Coord. of Hub Wall Pts.) (RHUB(I),LE1,NTIP) | 8F10.4 

IF NHUB>B = Gt (RHUB(1), 19, NTIP) | 8F10.4 
1F WHUB>16 = (G] (RHUGCIT), 117, NTIP) [| 8F10.4 
If WHUB>24 [G[ CRHUB(T),1=25, NTIP) {| 8F10.4 
IF MHUB>32  jG| (RHUB(T), 1=33,NTIP) | 8F10.4 
IF WLOSS>O, J=1 {G[ (Loss Parameter) (OLOS(K,J,1),K=1,5)  (Dfactor Array) (OFTABCK,J,1),K=1,5)] 10F8.4 
Jz2 |G] CDLOS(K, 5,1), K=1,5) (DFTAB(K,J,1),K=1,5)| 10FB.4 

- fel (DLOS(K,J,1),K=1,5) (DFTAB(K,J,1),K=1,5)] 10F8.4 

- fel (DLOS(K, 4,1), K=1,5) (DFTAB(K,J,1),K=1,5)| 10F8.4 

JENSTRM 1G] (DLOS(K,J,1),K=1,5) (DFTAB(K,J,1),K=2,5)| 10FB.4 

If MLOSS>1, J=1 |G] {Loss Parameter) (OLOS(K,J,2),K21,5)  (Ofactor Array) (OFTAB(K,J,2),K=1,5)| 10F8.4 
J=2 {6| (DLOS(K, 3,2), K=1,5) (COFTABCK,J,2),K=1,5)} 19F8.4 

- {cf COLOS(K,J,2),K21,5) (DFTAB(K,J,2),K=1,5)| 10F8.4 

- fa (OLOS(K,J,2),K=1,5) COFTAB(K,5,2),K=1,5)| 10FB.4 

JEWSTRK |G! COLOS(K,J,2),K#1,5) (DFTABCK,J,2),K=1,5)] 10°84 

IF MLOSS>2, J=t |G] (Loss Parameter) (DLOS(K,J,3),K=1,5)  (Ofactor Array) (DFTAB(K,J,3),K21,5)| 10F8.4 
J=2 |G| (DLOS(K,J,3),K21,5) (COFTABCK,J,3),K=1,5)] 10F8.4 

- fc] (DLOS(K, 5,3), K=1,5) (OFTAB(K,J,3),K=1,5)] 10F8.4 

- {S| COLOS(K,J,3),K#1,5) (DFTAB(K,3,3),K=1,5)} 10F8.4 

J=EMSTRM 1G] (OLOS(K,U,3),Ke1,5) (DFTAB(K,J,3),K=1,5)| 10F8.4 

IF MLOSS>3, J=1 [G] (Loss Parameter) (DLOS(K,J,4),K=1,5)  (Dfector Array) (OFTAB(K,J,4),K=1,5)] 10F8.4 
J=2 [G| (DLOS(K, 3,4), K=1,5) (COFTAB(K,3,4),Ke1,5)] 10F8.4 

- [sf (DLOS(K, J,4),K=1,5) (DFTABCK,J,4),Ke1,5)] 10F8.6 

- cl (DLOS(K,J,4),Ket,5) (DFTAB(K,J,49,K=%,5)] 10F8.4 

J=wsTam |G! COLOS(K,J3,4),K21,5) COFTABCK,J,4),K=1,5)] 10FB.4 

IF MLOSS>4, J=1 |G] (Loss Parameter) COLOS(K,J,5),K=1,5)  (Ofector Array) COFTABC(K,J,5),K=1,5)] 10°8.4 
J22 |G] (DLOS(K,J,5),€21,5) (DFTAB(K,J,5),K21,5)] 10F8.4 

- sf (DLOS(K,J,5),K=1,5) (OFTAB(K,J,5),K=1,5)] 10°8.4 

- {cj (DLOS(K,J,5),K=1,5) COFTABCK,J,5),K=1,5)|] 10FB8.4 

J=nsTRe [Gf (DLOS(K,J,5),K=1,5) (COFTAB(K,J,5),K=1,5)[ 10°8.4 


| Figure 20. General input data field format. 
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Sbtda cea ee Ra eee oS eee eS ee ee ee Se oS ee 2S. 6 8 6 Oe ee SS Ss eee ee Smee 


ieee Gi ee eee Se ee eS Bee eS Slates /6 ele Sree SS Re eR RE ST ae 


IF 


{7 | 


[yi RAR OvRe Mea Rs Resex cesar seN er oto tres Colum 


{R|ILOSS OPM OP OPO AA AB 
|RIBLADESCIRW) SOLIDCIRU) TILTCIRW) PRACERW) PRECIRU) PRCCIRW) PRDCIRW) 


Option jp| 
{al 

[A 

(a 

Jal 

{R | 

[R| 
ILOSS(1ROW)>=0{R | 
[R| 

|r| 

IDEF(IROW)#O = |R| 
IDEFCTROW)#O = [RI 
IDEFCTROW)#O {RI 
IDEFCIROW)#O = [R| 
IDEF(IROW)#O = [R| 
IDEF(TROW) MO = R] 
IOEFCIROW)#O = |R| 
IDEFCIROW)#O = {R| 
IDEFCIROW)#FO = |R| 
=TABLE jr] 
IF WSTRM>S = [R| 
BB=TABLE {rR | 
1F WSTPM>B = [R] 
CC=TABLE [®| 
IF WSTRM>B [R] 
DO=TAGLE {r| 
If WSTRM>B JR] 
EETABLE jr] 
IF NSTRM>8S [Rk] 
OPFVEL.DIA. = |R| 
OPeVEL.DIA. = {R] 


Input Data Field for 


the Program AXIDES 


(Data for Annular Stations & Rotors) 


SERS RETEST OR EHSEOA TASES OE ESSE SSS: 


VUTVAT119122222222223333 53333 366446644445 555555555 666666666 77777777778 | 
Je] 1234567890 123456 78901234567890 1234567890 1234567890 1234567890 12345678901 234567890 | 


AA ZT1P(1) ZHUB( 1) BI(1) 
AA ZTIPC(2) ZHUB( 2) BY{2) 
AA ZTIPC(3) ZRUB(3) BT(3) 
AA ZTIP(4) ZHUB(4) BI¢4) 
AA ZTIPC INAS) ZHUB(INAB) BIC1) 


DLIMCIRW) ALIMCTRW) BT(1) 


PTTCIRW) CPTC( J, 1RW), J21,5) 
TALECIRW) TBLECERW) TCLECIRW) TOLECIRW) 
TAMAX( ) TBMAX( ) TCMAX( ) TDMAX( ) 

(ACFE( 5S, TROW), J=1,4) 
(COR CJ, ROW), J=1,4) 
CACR(J, ROW), J=1,4) 
(CORI, ROW), J=1,4) 
(ELEC J, [ROM), J=1,4) 
CATECL, EROW), J=1,4) 
(CCTFCS, ROW), J=1,4) 
CATR I, DROW), J=1,4) 
CCTR(J, IROW), J=1,4) 
(incidence Angte Array on Streamlines) 


(Deviation Angle Array on Streamtines) 
Cintet/Outlet Segment Turning Rate) 
(Trans. Point Location) 

(Max. Thickness Point Location) 


(ZTEMPC(I-1,J),3=1,5) 
(ZTEMP(1,J),J3=1,5) 


BH(1) BLEEDC1) 
BH(2) BLEED(2) 
BH(3) BLEED(3) 
BH(4)} BLEEDC4) 
BH(I) BLEEDC!) 


cc DD EE 


BH(1) BLEED(I) CRENGYCIRBMATLCIR) NXCUTCIRW)| 
EB CHOKECIRW)| 
PRECIRW)| 


TATECIRW) TRTECIRW) TCTECIRW) TOTECIRW)| 
CHORDA( ) CHORDB( ) CHORDC( ) IDEF(IRW)| 


(BCF(S, 1ROW),J=1,4) 

(DCF (J, TROW), J=1,4) 
CBCR(J, ROW), J=1,4) 
(DCRCJ, 1ROW), 4=1,4) 
(ETECS, TROW), J=1,4) 
(BIFCS, (ROW), J=1,4) 
(DIFC J, 1ROW), J=1,4) 
(BIR(J, 1ROW), J=1,4) 
(DTR( J, IROW), J=1,4) 
CINCCTROW, J), J=1,NSTRM) 
CINCCIROW, J), J29,NSTRH) 
(DEV(IROW, J), J=1,NSTRM) 
(OEV(IROW, J), J=9,NSTRM) 
(PHI CTROW, J), d= 1, NSTRM) 
CPHICTROW, J), J=9, NSTRH) 
(TRANS (IROW, J), J=1,NSTRM) 
(TRANSC(IROW, J), J=9,NSTRM) 
(ZMAX(EROW, J), J=1,NSTRM) 
(2ZMAX(IROW, J), J=9, NSTRM) 
(RTEMP(1-1,8),021,5) 
(RTEMPC(I, J}, 3=1,5) 


ce sre ec ee re ee es ee 


ee re eh 


AS 6X, 5F10.4 
AG 6X, SF10.4 
AL, 6X, SFI0.4 
AS, 6X, 5F10.4 


AL, 6X, S5F10.4 
7F10.7,110 

15, 1X, 284,2x,3(AS, 
BF10.4 

6F10.4 

FIO. 
7F10. 
BFI0. 
BFIO. 
BF10. 
BF10. 
BF10. 
BF10. 
8F10. 
BFIO. 
BF10. 
BFIO. 


BF 10. 
8F10. 
8F10. 
8F10. 
BF 10. 
8F10. 
BF10. 
SFO. 


—_ 
oO 


- 


SB PF RP Rh PrP PF Be ER FAH AH BAA WA HH & & 


ee ee aa a Site iwiete la cae Saw waa e leew le ear Slo chee ewiao Se ee SoS 61 Se Slee ea wie C5 ale wel a See leis ww eee Sse he S SSS SS SSS ES SSeS Sm Se 


Figure 21. General input data field (cont'd). 
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Input Data Field for the Program AXIDES 


(Data for Stators & Annular Stations) 


IT] | 


Jy |**seeeneeceenccrecescnesesannvecnane Colum eee eG ene Cae e he sa enewes age ene) 


Option Ip] 191991117 122222222223333333333466466466455555555555666666665 77777777778 | Format 
| © | 1234567890123456789012345678901234567890 1234567890 1234567890 1234567890 1234567890 | 
|S AA ZTIP(INAB) ZHUB(INAB)  § BI(T) ss BHT) BLEED(I) | AG, 6X,5F10.4 
JS} OLIMCIRW) ALIMCTRYW) BT(1) BK(I) BLEED(1) NXCUTCIRW)| 7F10.7,110 
JS|#LOSS OPM OP OPO AA AB 88 cc DD £E EB CHOKECIRW)| 15,1K,2a4,2x,3(A4, 
|S|SLADESCIRW) SOLIDCERW) TILTCIRW) PRACIRW) PRECIRW) PRCCIRW) PROCIRW) PRECIRW)| 8F10.3 
1F ILOSSCEROW)>=0(S] PTTCIRW) (PIC(J, IRW), J=1,5) | 6F10.4 
JS] TALECIRW) TBLECIRW) TCLECIRW) TOLECIRW) TATEC(IRW) TBTECIRW) TCTECIRW) TOTECERW)| 8F10.4 
{S| TAMAX( ) TBMAX( ) TCMAX( ) TDMAX( ) CHORDA( ) CHORDB( ) CHOROC( ) IDEFC(IRW)| 7F10.4,110 
IF IDEF(IROW)#O S| CACF(S, ROW), J=1,4) (OCF(J, IROW), J=1,4) | BF10.6 
IF IDEFCIROW)#O S| (CCFC I, 2ROW), J=1,6) (DCF(S, ROW), J=1,4) | 8F10.6 
IF IDEFCIROW)#O |S} (ACR(J, FROM), J=1,4) (BCR(J, ROW), J=1,4) | 8F10.6 
1F IDEFCJROW)#O | S| (CCRC J, IROW) , J=1,4) (DCR J, ROW), J=1,4) | 8F10.6 
1f IDEFCIROW)MO |s| (ELEC J, IROW), J=1,4) (ETEC I, LROW), J=1,4) | 8F10.6 
If IDEFCIROW)#O {S| (ATFCJ, IROW), J1,4) (BTF(J, IROW), J=1,4) | 8F10.6 
if IDEFCIROW)#O [S| (CTFC(S, 1ROW), J=1,4) (OTFCS, IROW), J=1,4) j 8F10.6 
IF 1DEFCIROW)@O = |S{ CATRCJ, IROW) ,J=1,4) (BTR(J, IROW), 21,4) | 8F10.6 
1F 1OEFCIROW)MO |s| CCTRCS, FROW), J=1,4) (DTIR(J, IROW), Je, 4) | 8F10.6 
IF AA=TABLE |S] (Incidence Angle Array on Stream ines) CINCCIROW, J), J=1,NSTRM) | 8F10.4 
IF NSTRM>B S| CINCCTROM, J), J=9,NSTRM) | 8F10.4 
If @B=TABLE |S} (Deviation Angle Array on Streamlines) (OEV(TROW, J), J=1,NSTRM) {| 8F10.4 
IF WSTRM>B |S] (DEV(IROW, J), J=9, NSTRM) {| 8F10.4 
IF CC=TABLE |S] Cinlet/Ourttet Segment Turning Rate) (PHI CIROW, J), J=1,NSTRM) [ 8F10.4 
IF WSTRM>8 {S| CPHICIROW, J), J29,NSTRM) | 8F10.4 
IF DD=TABLE iS] (Trans. Point Location) (TRANS(IROW, J), J=1,NSTRM) | 8F10.4 
IF NSTRM>B [S| (TRANS (1ROW, J), J=9,NSTRM) | 8F10.4 
If EE=TABLE [S| (Max. Thickness Point Location) (2MAX(TROW, J), J=1,NSTRM) 1 8F10.4 
IF NSTRM>B S| (ZMAX(IROW, J), J29, NSTRM) | 8F10.4 
If OP2VEL.DIA. |S] (ZTEMP(I-1,2),J=1,5) (RTEMP(T-1,3),u=1,5) | 10F8.4 
Tf OPSVEL.OIA. = | S| (ZTEMP(1,J3), 3=1,5) CRTEMP(1, J), J¢1,5) | 108.4 


ELL ELE LL RAREST 8 SSS eS eS SO OOS OO eS SRN eH Sb Rew Kee e ewe ee ee eee eC iso ws wedewwesemecud cueeececdc Cceewsnwcececes 


atacand Te Oe 


[A] aa ZTIPCI) «= 2HUBCE) BTC I) BHC) BLEEDC!) | AS, 6x, 5F10.4 
Ja] aa ZTIPC1) © ZHUBCI) BTC1) BH(T) BLEEDCI) {| A4,6X,5F10.4 
fa] AA ZTIPCT) «© ZHUBC(T) BICt) BHC!) BLEEDCI) | A4,6x,5F10.4 
IF WXCUTCIROW)<O [C]  (Redial Location of Blade Sections) (XCUTCJ), J21,NC) | 8F10.4 
IF WC>8 |e] (XCUT(1), J®9,NC) | 8F10.4 
IF wC>16 = fc] (XCUT(J), JE17,NC) | BFI0.4 
{-| 
I-| 
|-] 
|-] 


Figure 22. General input data field (cont'd). 
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APPENDIX B. AXIDES OUTPUT 


This appendix contains the output file Examp.out, associated with input file 
Examp.dat. The input file is the example provided with the code when it was purchased 


and was the first to be successfully run on the AA Department workstations. 
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C60 °c 00% SZ- 000 SZ 


; 00% 'Sc- 
ZOE L 08h Ot- 00b'°Sz Ost Ot- 
Z0C TL 09S °St- 006 SZ 09S °SEt- 
fw) (uy) {Ww >) (w>) 
sntped 3 PUT PIOOD sntped e\eupprz0e) 
qoH petxy qnH dtl Teyxw dtL 
“SdUNOLNOD dNH GNW dId dOd SLINIOd WLY0 LOdNl 
e mre 7 ae Pang | T 
000°0 oO OCTOT GLT' 8Be Tl 
59900 °T G1 000°0 O° OLTOT OLT' 882 OT 
LSZ6 0 6 000°0 O'OCTOT OLT' 882 6 
Sts ‘6 e 000°0 O'oCTOT OLE ' 882 8 
G99.°0 u 000 °0 O'OLTOt OLT' 8982 L 
93899 '°6 9 000°0 O° OCTOT OLT' 882 9 
vits'o 4 000 °0 Oo OtTOl OLT' 882 S 
9690 0 4 000°0 O°OCTOr OLT' 882 b 
009 °6 t 000°0 O' Otto! OLT' 882 t 
Lobe 0 Zz 000°0 O'OCTOr OLT' 882 Z 
LSZ1°0 T 600 °0 O'OCTOl OLT' B82 T 
(99S/H) (W bs/6y) (°y °6eq) "ON 
ioTyseijg MOolg ae) g| AApOOTIA ainssalg ainqeiaduay aut lweardas 
syqnquesi4s aqnywee tis Tatuy yaTUI TEACL eT UY TeaoL Ie TUL 
AGOLIWBSULS YO ANITWWAYLS AG SNOILNSIYLSIC Lndil 
© 
O 
teebe SI 3CbbS7Z°0 ¢ PeodeZ7i AbG6TB G- * Coede60-S31S918°O + ZeodeID-FLIPBS O- * LePO-3LSS6£'°0 + OOTA7T9ILEZ O > d> 
w1d4 Burmoy, [Og AuL ut st petwoudtog 3esaH IFJTIedS SUL 
-suotqyeys rertnuuy g ‘3 PUY sa6pyZ apetgd euL IF ePeEW eG TI tM SUOFABTNOTeO 
“smaoy epeTad 2 SPH rosseidwoy auL ‘sauy[wesays T] uo pawzojieg eq [TIN SuofFAe[NoTeod 
6°82 St AHTay TeTNOaTOH aUL ' p2S't Ss} OF ey aeansseid xossailduoD paeijsed euyuL 
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APPENDIX C. TEST RIG 


This appendix contains information pertaining to the transonic axial compressor 
test rig located at the NPS TPL. Included is the general procedure involved in operating 
the test rig, and its disassembly and reassembly. This information is provided to fill the 
void that currently exists in the documentation associated with maintenance and operation 


of the test rig. 


A. OPERATING PROCEDURE 


A photograph of the control panel, from which the operation of the test rig is 
controlled and monitored, is located in Figure 23. The panel is used to control the supply 
of air from the laboratory’s 12 stage compressor to the test-rig drive turbine, to relieve the 
axial load on the drive shaft bearings by porting low pressure air to a balance piston inside 
the compressor hub, and to remotely operate the hydraulic throttle located in the 
compressor-intake housing. It is also used to monitor the Operating temperatures of the 
drive-shaft bearings, the RPM, and levels of vibrations. The following is a description of 


the operating procedure. 


1. The technician starts and warms up the laboratory air-supply compressor (1 
hour), with the bypass (dump) valves in the open position. | 

2. Turn on the control-panel activation switch (located beneath the control-panel 
counter-top). 

3. Ensure that air is being supplied to the balance piston, and bearing oil-mist 
cooling lines. This air is supplied through the shop air lines, by the Elliot 
compressor. 

4. Set the oil-mist cooling, air-oil regulators, to 35 psi, and 12 drops of oil per 
minute (count the number of drips over a three minute interval and average). 

5. Check that the bearing temperature readings are indicating test-cell ambient 


values (not wandering). 


93 


6. Close supply-air dump valves (2 ea.),until a two-to-one pressure ratio across 
the laboratory air supply compressor 1s indicated (on gauge to the right of 
control panel), and the noise level in the control room is acceptable (some 
combinations of valve settings are noisy). 

7. Open air supply valve and close dump valves partially until desired RPM is 
indicated, while maintaining the two-to-one pressure ratio across the air 
supply compressor. 

8 As the RPM increases, the air supplied to the balance piston must be increased 
to compensate for the axial load on the bearings. If the value displayed on the 
analogue readout is maintained at the level indicated before start-up (with no 
air supplied to the balance piston), the axial load produced by the spinning 
rotor will be negated by the pressure on the balance piston. This will result in 
approximately zero axial load on the high speed bearings. 


9. Set throttle valve as required, and correct air-supply valve to maintain RPM. 


During operation of the test rig, RPM, bearing temperature, balance piston, and 
vibration indications must be monitored continuously. Additionally, periodic checks of the 


oil-mist regulators is recommended. 


B. DISASSEMBLY AND REASSEMBLY 


The disassembly and reassembly performed in association with this report was 
documented on video tape. This tape will be maintained with the records of the test rig for 
guidance in future overhaul efforts. A few notes, however, concerning the disassembly and 


reassembly of the machine should be emphasized. 
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Disassembly 


The removal of the inlet nozzle bell from the hub housing is tricky. Extreme care must 
be taken to prevent the shroud (attached to, and removed with, the nozzle bell), from 
coming into contact with the rotor blade tips. To remove the nozzle bell and shroud 
assembly without incident, pull out the positioning pins, and un-screw the bolts 
securing the assembly to its support structure. Then attach the test-cell overhead crane 
to the eye bolt on the top of the nozzle and raise the crane until there is a slight tension 
on the chain. Attach a winch to the nozzle-inlet face, and the piping from the inlet 
housing (outside the test cell), and gently pull it axially away from the hub housing. If 
the winch cables are evenly distributed around the nozzle face the assembly should 
slide along the channels fastened to the nozzle-bell support structure that ensure 
clearance of the rotor blade tips. The overhead crane will prevent the assembly from 


dropping as it clears the edge of the nozzle base. 


Reassembly 


Ensure that markings on splines are located during disassembly and aligned during 
reassembly. These markings line up the high speed rotational components as they were 
when the machine was originally balanced. 

The pre-load on the high speed bearings is defined as “hight”. A mere 2-3 in.-lb. 
torque, on the rotor mounting bolt, is all that is desired to keep the bearings, rotor and 
rotor shaft in place. 

Perform a new calibration of the torque balance each time the machine is reassembled. 
Check bearing  oil-mist cooling lines, bearing temperatures, and _ internal 


instrumentation lines for proper operation/indications prior to reassembling. 
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C. INSTRUMENTATION 


The test runs and overhaul of the test rig, as described in this report, have been 
primarily driven by the desire to return the machine to routine operation. The myriad of 
instrumentation, that was added in a piece-meal fashion over the fifteen years of operation 
of the rig, is in need of overhaul. Currently the data-acquisition system at the TPL is in 
transition. The Hewlett-Packard Basic workstation and programs that have been used for 
data acquisition are being augmented and possibly replaced with a PC data acquisition 
controller using Labview™. A VXI-Bus mainframe has been purchased, and obsolete 
scanners are to be replaced. All instrumentation lines from the test-cell to the data 
acquisition area must yet be verified or replaced before performance mapping of the new 


stage can begin. 
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Figure 23. Transonic Compressor Test Rig Control Panel. 
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